Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 



(MN M?eX style file v2.2) 



Towards a Resolution of the Galactic Spin Crisis: 
Mergers, Feedback, and Spin Segregation 

Ariyeh H. Mailer and Avishai Dekel 

Racah Institute for Physics, The Hebrew University, Jerusalem 91904, Israel 



1 February 2008 



ABSTRACT 

We model in simple terms the angular-momentum problems of galaxy formation in 
CDM cosmologies, and identify the key elements of a scenario that may solve them. 
The buildup of angular momentum is modeled via dynamical friction and tidal strip- 
ping in a sequence of mergers. We demonstrate how over-cooling in incoming halos 
leads to a transfer of angular momentum from the baryons to the dark matter, in con- 
flict with observations. By incorporating a simple recipe of supernova feedback, we are 
able to solve the problems of angular momentum in disk formation. Gas removal from 
the numerous small incoming halos, which merge to become the low specific angular 
momentum (j) component of the product, eliminates the low-j baryons. Heating and 
puffing-up of the gas in larger incoming halos, combined with efficient tidal stripping, 
reduces the angular momentum loss of baryons due to dynamical friction. Dependence 
of the feedback effects on the progenitor halo mass implies that the spin of baryons 
is typically higher for lower mass halos. The observed low baryonic fraction in dwarf 
galaxies is used to calibrate the characteristic velocity associated with supernova feed- 
back, yielding Vfb ~ 100 km s""'^, within the range of theoretical expectations. We then 
find that the model naturally produces the observed distribution of the spin parameter 
among dwarf and bright disk galaxies, as well as the j profile inside these galaxies. 
This suggests that the model indeed captures the main features of a full scenario for 
resolving the spin crisis. 
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1 INTRODUCTION 



The 'standard' model of cosmology, which assumes hierar- 
chical I 



: uildup of structure in a universe where the mass is 
dominated by cold datk iiiaLLet (CDM), seems Lu be fa 



tac- 



ing intriguing difficulties in explaining some of the robust 
observed properties of galaxies. Standing out among these 
problems is the inability of galaxy formation models to re- 
produce both the sizes and structure of disk galaxies. In hy- 



drodyn^mical simulations baryons lose a signmcant traction 
of their atlgukr ttiotilcitltrltil leading to ovCrly sfllall diriks. 
To avoid this, analytic and semi-analytic models commonly 
assume there is no angular momentum loss. Recent studies, 
however, show that this leads to a variety of discrepancies 
with observations. Taken together this suggests a crisis in 
our understanding of the role of angular momentum in form- 
ing disk galaxies. Our aim here is to make progress in the 
effort to resolve this crisis by first reproducing the problems 
using a simple model in which the important physical in- 
gredien ts are spelled out and well understood. Subsequently 
we incorporate an additional process into this model which 



then reproduces the observed sizes and structure of galactic 
disks. 

The sizes of galactic disks are commonly linked to the 
angular moment um of their parent dark-matter halos (Fall & 



Efstathiou 1980). This modeling is based on the distribution 



of halo spin parame ters as found in N-body simulations (see 
Bullock et al. 2001, and references therein), combined with 



the assumptions that the baryons and dark matter initially 
share the same distribution of specific angul ar momentum. 
7, within th e halos (as seen in simulations by van den Bosch 



et al. 2002) and that j is conserv ed as the bar yons contract 



to form the disk (as suggested by Mestel 1963). The sizes of 
disks obtained under these assumptions are roughly compa- 
rable to those observed. 

However, high-resolution simulations that incorporate 
gas processes find this scenario to be invalid. In particular, 
they find that a significant fraction of the angular momen- 
tum of the baryons is transfered to the dark matter, result- 
ing in disk sizes roughly an order of magnitude smaller then 



et al. 



those observed ( Navarro fc Steinmetz 200C jommer-Larsen 
I999J [Navarro fc Steinmetz 1997| ; [Navarro et al. 199. 
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Navarro & Benz 1991). This has been refered to as the an- 



higher spin than average (Gardner 2001; Wechsler 2001) 



gular momentum catastrophe. 

The angular momentum catastrophe is commonly asso- 
ciated with the problem of "oyer-cooling" in CDM -type sce- 
narios ([White fc Frenk 1991^ [White fc Rees 1978|). Without 



sufficient feedback much of the gas cools quickly, contracts 
into small halos and then spirals deep into the centers of 
bigger halos, efficiently transf erring its orbital angular mo - 
mentum to the dark matter (Navarro & Steinmetz 2000). 



It has therefore been speculated that enough energ y feed- 



bac k from supernova may p revent this over-cooling (Larson 
1974; \^^hite & Rees 1978) and thus reduce the angular- 



momentum loss. Indeed, simulations where gas cooling is ar- 
tificially suppressed till z = 1 do not suffer from the angula r 
momentum catastrophe (Eke et al. 2000; Weil et al. 1998). 
Furthermore, simulations including some forms of feedback 
s how a reduction in the amo unt of angular momentum lost 
( 3ommer-Larsen et al. 1999). However, while feed back has 
been studied using simplified approximations (e.g., Dekel & 
Silk W M. hereafter DSI, 



back ha s proved challenging (see ^^^^^^^^^^^^^^^^^ 

and references therein), though some partial progress may 
have been made r ecently ( Ferrara & Tolstoy 2000 ; Thacker & 
Couchrian 2001; Springe! fc Hernquist 20Q2[ ). At this point. 



the feedback scenario has not yet been studied in satisfactory 
detail, nor has it been confirmed to solve the spin problem, 
or properly understood in basic terms. This motivates an at- 
tempt to understand the scenario and how it may work using 
a simple semi-analytic model, in which the basic elements 
are easily understood. 

Even if the angular momentum catastrophe can be 
avoided such that the total baryonic spin agrees with the 
total dark-halo spin, there remain discrepancies between the 
angular-momentum properties of dark halos in N-body simu- 
lations and those of observed disk galaxies. Foremost among 
these is the apparent mismatch of the distribution of specific 
angular momentum within galaxies (loosely termed "t he j 
profile" ) between observations and N-body simulations. Bi 
lock et al. (2001, hereafter BD) revealed a universal pro- 



file with an excess of both low-j and high-j material com- 
pared to the observed j distribution in exponential galac- 



tic disks. Subsequently, van den Bosch, Burkert & Swater; 



(200lQiereafter BBS) have measured the j distribution in 
a sample of dwarf disk galaxies, confirming in detail, case 
by case, the differences between the j profile predicted by 
the simulated halos and those found in disk galaxies. We 
refer to this second discrepancy as the mismatch of angular- 
momentum profiles. 

Another spin problem is indicated by the observations 
of de Jong & Lacey (2000), who found that the scatter in 



the log of the spin parameter inferred from fitting the spread 
of observed disks sizes in a large sample of late-type disk 
galaxies is only 0.36 ± 0.03, while the scatter seen in sim- 
ulated halos is significantly larger, 0.5 ± 0.05 (BD and ref- 
erences therein). This could have been explained by a sce- 
nario in which the disks formed in low-spin halos are un- 
stable and thus become early-type galaxies, which implies 



that di sk galaxies occupy only the high er spin halos (Mo 



Mao & [White 1998^ [van den Bosch 1998)) . However, this sce- 
nario is in apparent confiict with the finding in N-body sim- 
ulations that recent major mergers, usually identified with 
large spheroidal components, actually give rise to halos with 



This demonstrates that the exact connection between the 
angular-momentum properties observed in galaxies and in 
N-body simulations is unclear. 

We propose that the solution to the crisis is spin seg- 
regation — that the angular momentum distribution of 
baryons differs from that of the dark matter due to gas pro- 
cesses. These processes can either decrease or increase the 
specific angular momentum of the baryons relative to the 
dark matter. Gas cooling generally results in lower spin for 
the baryons compared to the dark matter, but heating due 
to feedback reduces this effect, and gas removal in small ha- 
los actually results in a higher spin for the baryons compared 
to the dark matter. 

In this paper we work out a simple model to explore 
these effects. Knowing that in a hierarchical scenario the 
halo buildup can be largely interpreted as a sequence of 
mergers between smaller halos, our model is based on a sim- 
ple algorithm for the buildup of halo spin by summ ation of 
the orbital angular momenta of merging satellites (Mailer 



et al. 2002; Vitvitska et al. 2001). This algorithm has been 



found to match well the spin distribution among halos in 
N-body simulations. In this paper we extend the model and 
find that it also reproduces the angular-momentum profile 
within halos. It therefore provides a useful insight into the 
origin of the spin problems, and a clue for their possible solu- 
tion. We therefore use this model as a tool for incorporating 
the relevant baryonic processes, and especially feedback. 

In ^ we describe the model for angular-momentum 
buildup in halos by mergers. In §^ we model the angular 
momentum catastrophe as resulting from over-cooling and 
dynamical friction. In ^ we introduce a simple model of 
feedback. In §^ we study the resultant angular-momentum 
properties of the baryonic component both in bright and 
dwarf disk galaxies and compare them to observations. In 
^ we test the robustness of our results to the details of 
the feedback model assumed. We conclude and discuss our 
results in 



2 BUILDUP OF HALO SPIN BY MERGERS 

2.1 Modeling the Spin Parameter 

The angular momentum of a galaxy, J, is commonly ex- 
pressed in terms of the dimensionless spin parame ter A — 
JyJjE\/GM^^^ , where E is the internal energy (Peebles 
196£ ). In practice, the computation and measurement of this 
quantity, especially the energy, may be ambiguous, e.g., it is 
not obvious how to define separate spin parameters for the 
dark matter and the baryons. Furthermore, it introduces an 
undesired dependence on the specifics of the halo density 
profile. Instead, following BD, we use the modified spin pa- 
rameter 



A' = 



\/2Vvir^vir 



(1) 



where j = ,J/M is the specific angular momentum. This 
quantity is more straightforward to compute for each com- 
ponent, \'^^ and A{,, and it does not explicitly depend on 
the density profile of the halo. For a d ensitv profile of the 
type suggested by Navarro et al. (1996 , NFW) A' is related 
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to A by 

A' = Xf-'^' . 



(2) 



where the factor fc reflects the difference between the energy 
of an NFW halo and a truncated singular isothermal sphere. 
An expression f or fc as a functio n of the halo concentration, 



Mo et al. (1998, equation 22). For a NFW 



4.5 the values of the two spin parameters 



Cvir, IS given m 
halo with Cvir - 
become equal. 

BD found for simulated halos that, similarly to A, the 
distribution of A' is well fit by a log-normal function, 



P(A')dA' 



1 



^2^ 



exp 



2al 



dX' 



(3) 



with A n 0.035 (compare d to Ao ~ 0.042) and cta ^ 0.5. 



In 



Mailer et al. (2002), we proposed a very simple model 

the orbital- 



for the buildup of spin in dark halos, termed 
merger model" . In this model, the spin of a flnal halo is sim- 
ply the vector sum of the orbital angular momenta of all the 
halos merging into its main progenitor th roughout its his- 



tory. In this paper, as in Mailer et al. (2002 ) , we generate 500 



random realizations of merger histories for each desired flnal 
halo mass at z = 0, using the method o f [Somcrvillc fc Kolatt 
(1999^ with the slight adjustment of Bullock, Kravtsov & 
Weinberg ( 200C ) to implement the Extended Press Schechter 
formalism ( [Lacey fc Cole 1993 Bond et al. 1991 ). The spe- 



ciflc cosmology assumed, without significant loss of gener- 
ality, is the "standard" ACDM cosmology, with fim = 0.3, 
Q.A = 0.7, h — 0.7 and as = 1.0. For each of the merger trees, 
we perform 10 realizations of the orbital angular momentum 
that comes in with each merger. The encounter parameters 
are taken to be those of typical encounters in simulations, 
and the directions of the added spins are drawn at random, 
allowing for a slight correlation between t he planes of succes- 



sive mergers (as see n in simulations, e.g. Dekel et al. 2001 



Porciani et al. 2001). We end up with a distribution of 5000 
spin values for each halo mass. This distribution is shown 
in Fig. ^ in comparison with the log-normal distribution 
obtained by BD from a cosmological simulation. The model 
predictions of Aq = 0.036 and a\ = 0.53 matc h rather well 



the simulated results. One can see in Figure 2 of Mailer et al 



(2002^hat the match is reasonable even without introduc- 
ing any correlation between the enco unter planes 



Similar results were obtained by Vitvitska et al. (2001) 



who draw their encounter parameters statistically from the 



distribution found in N-body simulations, and by Monaco 



Theuns fc Taffoni (2002), where encounter parameters for 



each object are provided by Lagrangian perturbative the- 
ory. We conclude that the simple orbital-merger model is 
successful in reproducing the simulated halo spin distribu- 
tion, and therefore serves as a useful tool for studying the 
angular-momentum buildup in galaxies. 



2.2 Modeling the Angular-Momentum Profile 

The cumulative mass distribution of specific angular mo- 
mentum within simulated dark halos, M{< j), was found 
by BD to be well fit by the 2-parameter universal functional 
form 



M{< j) 



M - 1 +j/jn 



At > 1, 



(4) 
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Figure 1. Dark-halo angular-momentum profiles: The mass frac- 
tion with specific angular momentum less then j as a function of 
i/imax. The lines connecting symbols represent a random sample 
of profiles produced by our orbital-merger model. The smooth 
curves are fits to these profiles, where /i is chosen to cross at 
half the total mass. The shaded region marks the 90% spread of 
profiles found by BD in simulations. 



in which one of the parameters, say j'max, can be replaced 
by A'. This profile is a simple power law, M{< j) oc j, for 
at least half the mass. When fi is close to unity, it bends 
over towards the high-j end. The distribution of the shape 
parameter /i was found to be Gaussian in log (/-t — 1), with 
a mean of —0.6 and standard deviation of 0.4. 

One can also quantify the shape of the angular momen- 
tum profile by the ratio of the total specific angular momen- 
tum, jtot, to jmax. The relation between this parameter and 
fi is 



Jtot 
Jmax 



Kl-(l-(ft-l)log 



fI-1 



)) 



(5) 



(BBS eqn. 12). 

We create M{< j) profiles for each of the realizations 
of the orbital-merger model as follows. We divide the mass 
growth of the halo into 20 equal-mass bins and assign to 
each bin the corresponding angular momentum that comes 
in with that mass. When a satellite's mass is divided among 
several bins, i = l,n, with a fraction / in bin 1, we assign a 
fraction of its angular momentum to bin 1. The remaining 
satellite mass and angular momentum are then distributed 
in an analogous way among bins i — 2,n. This mimics an 



M{< j) 



oc J 



1/2 



profile for the contribution of each satel- 



lite, in rough agreement with the predictions of a toy model 
based on dy namical friction an d tidal stripping for an NFW 
halo (BD; Dekel et al. 2001). The role of this procedure 



is to smooth the profile; in most cases it does not change 
the global shape, which is predominantly determined by the 
cosmological sequence of mergers. The result is therefore not 
too sensitive to the exact assumed distribution of j due to 
each single merger. 

Fig. hi shows a sample of j profiles produced by this pro- 
cedure. Our profiles are somewhat noisier then those seen in 
simulations, partly because we ignore any exchange of an- 
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Figure 2. The distribution of sliape parameter ^ for halo angular- 



momentum prohlcs. i'he histogram describes the distribution in 
the sample reproduced by our orbital-merger model. The smooth 
curve is the log-normal fit to the N-body simulation results of 
BD. There is a slight offset in the mean of the distribution of our 
model (-0.5 instead of -0.6) but this is not statistically significant. 



gular momentum between particles. The overall agreement 
between our model profiles and the simulation results of BD 
suggests that the rearrangement of angular momentum only 
smoothes the profile without changing its global shape. We 
fit each of our profiles to eq. by the simple method of re- 
quiring an exact match at half the virial mass (Fig. I). The 
same results follow if we use to compute /j. by inverting 
eqn. ^. We obtain a distribution of u, values that are consis- 
tent with the results of BD (Fig. Q). As in BD, they show 
no dependence on halo mass or redshift. Thus, the orbital- 
merger model, which reproduces the correct distribution of 
spin values among halos, also recovers reasonably well the 
angular-momentum profile inside halos. 

Our orbital-merger model reveals an important feature 
in the buildup of halo angular momentum. We see in the 
typical realization that the final halo spin and its direction 
are predominantly determined by one last major merger. 
Most of the high-j mass in the final distribution comes in 
with this merger. The many smaller satellites, each possibly 
carrying a high j, come in at different directions and there- 
fore tend to sum up to a low spin. Thus, most of the low-j 
fraction of the halo mass originates in minor mergers. This 
feature provides an important clue for a possible solution 
to the spin crisis, as follows. If small satellites lose a large 
fraction of their gas before they merge into the halo then 
the final gas distribution in the galaxy lacks the low j ma- 
terial. Furthermore, if more of the galactic gas originates in 
big satellites then the spin of the baryonic component may 
be higher then that of the dark matter. 



REPRODUCING THE ANGULAR 
MOMENTUM CATASTROPHE 



3.1 Tidal Stripping and Dynamical Friction 

We can understand the transfer of angular momentum from 
the baryons to the dark matter in cosmological hydrodynam- 
ical simulations with a simple model including gas cooling, 
dynamical friction and tidal stripping. In our model for halo 
spin buildup discussed above we have not specified how the 
orbital angular momentum of the satellite is converted to 
the spin of the halo. Clearly the transfer of angular momen- 
tum occurs by the processes of dynamical friction and tidal 
stripping. The dynamical friction acting by the halo parti- 
cles on the bound part of the satellite exerts a torque which 
transfers angular-momentum from the satellite to the halo 
and eventually brings the satellite towards the halo center. 
Satellite particles that are tidally stripped and become part 
of the halo before the satellite sinks to the halo's center re- 
tain what is left of their angular momentum at that point 
and add it directly to the halo. 

A hint for the typical M(< j) distributi on due to a 
singl e merger can be obtained as follows (BD; Dekel et al 



2001). If the satellite moves with a circular velocity Vc{r), 
and we assume that the lost mass and j are deposited, on 
average, locally at r, then the resultant spatial profile j{r) 
can be obtained by averaging over shells. 



4Tvr p(r)j(r) = m(r) 



+ 



dm(r) 
dr 



rVc(r) 



(6) 



The first term refers to the J transfer as a reaction to dynam- 
ical friction, and the second is the J deposit associated with 
the tidal mass transfer. We need to estimate the momen- 
tary surviving satellite mass m(r) in order to compute j(r), 
which can then be inverted to obtain the desired M{< j). 

The tidal mass loss at halo radius r can be crudely 
estimated by evaluating the tidal radius It of the satellite at 
r using the crude resonance condition. 



(7) 



where m{£) and M(r) are t he mass profile s of the satellite 
and halo respectively (e. g. Weinberg 1994). If, for simplic- 
ity, these two profiles are self-similar, then the resonance 
condition implies 



et/e.ir = r/R, 



(8) 



where ^vir and -Rvir are the virial radii of the satellite and 
the halo respectively. This means that the remaining bound 
mass of the satellite when it is at r obeys 



m[£t{r)] oc M{r). 



(9) 



A more accurate recipe for tidal stripping, based on studies 
with N-body simulations of mergers, reveal that eq. is 
valid as a crude approximation for a wide range of merger 
parameters (Dekel et al., in preparation). We therefore adopt 
it in our model. 

If the halo density profile is isothermal, M(r) oc r, then 
substituting eq. (^ in eq. (^) yields j{r) oc r, which corre- 
sponds to M(< j) oc j. A more realistic NFW density profile 
leads in the outer halo t o a j profile closer to M(< j) oc j^''^, 
(BD; Dekel et al. 2001) which is what we use to smooth the 
profile in the orbital-merger algorithm. The qualitative sim- 
ilarity between the profile predicted by this toy model and 
the typical j profile seen in simulations (BD) provides fur- 
ther support to the general picture of J buildup by mergers. 
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3.2 Spin Segregation 

We can now explore the effect of cooling in the satellites 
on the final angular momentum of the baryons. We assume 
that at an early time the baryons follow the dark-matter 
distribution. As the gas within a halo cools radiatively, it 
contracts to a more compact configuration characterized by 
a smaller radius, Rh, compared to the extent of the dark 
matter, iidm- Such a spatial segregation in a satellite before 
it merges with a bigger halo naturally leads to a difference 
in the spins of baryons and dark matter in the product halo. 
The bound part of the satellite is continuously transferring 
orbital angular momentum to the dark matter by dynamical 
friction. In the extreme case, the inner satellite that remains 
bound all the way to the halo center eventually loses all 
its orbital angular momentum by this mechanism. At the 
moment of stripping, the escaping satellite material stops 
being affected by dynamical friction, and it joins the halo 
with the angular momentum that it is left with at that point. 

The spatial segregation in the satellite thus implies that 
the j-rich mass stripped at the early stages of the merger in 
the outer halo is dominated by dark matter, while the more 
compact baryonic component survives longer as a bound 
satellite and loses more of its orbital angular momentum 
by dynamical friction. The result is a net transfer from the 
baryons to the dark-matter component. This process is il- 
lustrated in a schematic diagram. Fig. ^ 

Using eqn. ^ above, we obtain that the orbital angular 
momentum added to the merger product with the baryons 
is related to the amount gained by the dark matter via 

AJb = ^AJdn.. (10) 
Jtdm 

We ignore any transfer of angular momentum from the satel- 
lite to the baryons already in the halo as baryons in the more 
massive halo are already more centrally concentrated than 
those in the satellite. In the case of maximum cooling the 
baryons dominate the center of the halo which we approx- 
imate as Rh = fhRdm, whe re fb is the unive rsal baryon 
fraction. Adopting /b — 0.13 (Tytler et al. 1999), we obtain 
that the spin of the baryons is reduced by almost an order 
of magnitude, and thus reproduces the angular momentum 
catastrophe as seen in hydrodynamical cosmological simula- 
tions. In Fig. ^ we show the resultant probability distribu- 
tion of the spin parameter as obtained in our orbital-merger 
model with maximum cooling; it has Ag — 0.005. 

This simple modeling of the spin problem makes it easy 
to understand why it is usually assumed that some form of 
energy input into the gas may remedy the problem; it would 
delay the cooling, increase i?b, and thus reduce the baryonic 
spin loss. 

Note that in common semi-analytic models of galaxy 
formation the gas is assumed to be either in a hot phase or 
in a cold phase, with the hot phase extending out to -Rdm 
and the cold phase condensed. This implies that the cold 
gas in the satellite looses most of its angular momentum 
to dynamical friction, while the hot gas retains most of its 
angular momentum. In this case the fraction of baryons in 
the hot phase, fhot, determines the orbital angular momen- 
tum added as baryonic angular momentum to the merger 
product, AJb = /hotAJdm- Thus, one can relate our dis- 
cussion here to other semi-analytic modeling by identifying 

fhot = Rh/Rdm- 



Over-Cooling — > Spin Crisis 



Gas Cooling 




Figure 3. A schematic illustration of how over-cooling in merg- 
ing satellites leads to the angular momentum catastrophe. The gas 
contraction within the incoming satellite makes the gas immune 
against tidal stripping: it spirals all the way into the halo center 
while losing all of its orbital angular momentum to the dark halo 
due to dynamical friction. The dark matter, which dominates the 
outer regions of the satellite, is gradually stripped in the outer 
parts of the halo, thus retaining part of its orbital angular mo- 
mentum. 



4 FEEDBACK 

It is generally assumed that some form of heating may play 
an important role in preventing over-cooling and therefore 
angular-momentum loss. The most important source of heat- 
ing is likely to be supernova feedback though reionization by 
the UV background, tidal heating, and ram pressure may all 
contribute. Semi-analytic models assume different ad hoc 
recipes for the amount of heating as a function of star- 
formation rate, though hydrodynamical simulations have 
been largely unsuccessful so far in i mplanting feedback in 



a way that prevents over-cooling (see Thacker & Couchman 



200C, and references therein). Our approach in the present 



work is to avoid the details of star formation and feedback 
recipes and rather use a very simple prescription for the ef- 
fect of feedback in a satellite halo as a function of its virial 
velocity, Kat. 

Following the analysis in §IV of DS, we assume that the 
total supernova energy pumped into the gas per unit mass is 
a robust quantity. This is based on the assumptions that the 
supernova rate is roughly proportional to the same gas mass 
that is later being affected by it, that the star-formation rate 
is inversely proportional to the dynamical time, and that 
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Figure 4. The effect of over-cooling on tlie spin distribution 
of baryons compared to the dark matter. Since the baryons are 
tightly bound in the satellite centers, they spiral into the inner 
part of the big halo without being stripped and thus loose most of 
their angular momentum to the halo by dynamical friction. The 
model reproduces the deficiency of angular momentum in baryons 
of almost an order of magnitude as seen in hydrodynamical sim- 
ulations. 

the supernova remnant transfers most of its energy to the 
interstellar medium by the end of its "adiabatic" phase. DS 
estimated as a crude upper limit that the supernova input is 
comparable to the kinetic energy per unit mass of an isother- 
mal halo with virial velocity ~ lOOkms"^. The actual free 
parameter characterizing our feedback model, Vib, is defined 
as the virial velocity of a halo for which the energy input is 
sufficient to heat all the gas to the virial temperature of the 
halo. Based on DS we expect Vfh to be of order lOOkms"^ 
or less, but we keep this estimate in mind only as a general 
constraint which does not enter our current analysis. The ac- 
tual relation between Vfb and the supernova input depends 
on the initial state of the gas before the supernova burst. 

We thus parameterize the effect of feedback on the spa- 
tial extent of the baryons in a halo by the ratio of Vaat to 
Vfh- The limit Vsat >> V'fb, massive halos with deep poten- 
tial wells, corresponds to maximum cooling, 7?b <C ^dm- For 
smaller halos with I4at — Vfh we expect the heating to bal- 
ance the cooling and yield -Rb — Rdm- Our model is therefore 
a simple interpolation between these limits, 

Rh = (:^) Rdn. (Kat > T4b) , (11) 

with 71 an arbitrary exponent. 

If Vfb is larger than Vsat, the feedback can cause gas 
ejection from the satellite. For I4at << Vfb we expect total 
blowout, and we assume that partial blowout starts occur- 
ring for halos where the potential energy of the gas is compa- 
rable to the energy input from supernova, Vg^t = ( 1/2) Vfb- 
We therefore parameterize the amount of gas that remains 
in the halo by another simple interpolation, 

/d = {^^^y^ (^-t < Vfb/V2) , (12) 



0.15 




Vsatt (km/s) 

Figure 5. The effect of feedback on the ratio of baryonic angular 
momentum to that of the dark matter as a function of the virial 
velocity of the merging satellite. In the model described by the 
thick curve, the feedback is characterized by Vfb = 95kms~^, 
with the exponents 71 = 72 = 1- If the baryons in the satel- 
lite are distributed like the dark matter (middle range), then 
Jh/Jdm = /bi the cosmological baryon fraction. In larger satel- 
lites with Vsat > Vfb, the baryons are less extended then the dark 
matter and thus loose more angular momentum by dynamical 
friction before they are stripped into the halo. In smaller satel- 
lites with Viat < Vfb/V2i gas is blown out and therefore Jb is 
reduced, though the specific angular momentum of the baryons 
remains the same as that of the dark matter. The shaded region 
shows the range of models we consider in § ^ all constrained to 
match the baryonic fraction in dwarf galaxies. 



with 72 another arbitrary exponent. 

We start with the simplest arbitrary choice 71 = 72 = 1- 
Later, in section ^, we will explore the robustness of the 
results to different assumptions about choices of values for 
71 and 72. 

The only remaining free parameter of the feedback 
model is Vb- We set it by requiring that halos with a virial 
velocity of 60kms~^ have an average /d = 0.04, as mea- 
sured in the data of BBS. Note that we are concerned with 
the state of the baryons only when the satellite mergers with 
the halo; it is therefore sufficient that the feedback be effec- 
tive for a short duration, perhaps triggered by the merging 
process itself. 

The effects of heating and blowout on the angular mo- 
mentum of baryons versus dark matter, according to our 
model, in shown in Figure |^. For massive satellite with 
Vsat > Vfb, the baryons retain less angular momentum as the 
mass increases, because of dynamical friction. Small satel- 
lites with Vat < V^Vfb allow a blowout of a larger fraction 
of their gas as the mass decreases, thus adding to the merger 
product less baryons and therefore less baryonic angular mo- 
mentum, while the specific angular momentum of the two 
components remains the same. 

In Fig. ^ we demonstrate the effects of the different in- 
gredients of this feedback scheme, with Vfb ~ 95kms~^, on 
the baryonic spin-parameter distribution. We do it for two 
kinds of final halo masses, corresponding to virial velocities 
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Figure 6. The effects of heating and blowout in our model on the distribution of baryonic spin parameter A'. Log-normal fits are shown 
for each histogram, with the corresponding mean and scatter quoted. In each panel the higher peaked curve and higher two numbers are 
the bright galaxies (Vvir = 220kms~^) and the lower peaked curve and lower two numbers are the dwarf galaxies (Vvir = 60kms~^). 
Left: Only heating is included in our feedback recipe (blowout ignored). Center: Only blowout is included (cooling ignored). Right: 
The whole feedback recipe is included, with cooling, heating and blowout. 



of 60 and 220 km . We refer to these as representing dwarf 
and bright galaxies respectively. 

The left panel of Fig. |^ shows the effects of adding 
heating only, ignoring blowout. The dwarf galaxies are all 
built up by satellites of Vsat < V^h, so there is full heat- 
ing, 7?b — Rdm The result is that the spin distribution of 
the baryons resembles that of the dark-matter. In bright 
galaxies, the partial heating raises the baryonic spin to 
Ag = 0.028, which is still lower than the typical spin of the 
dark matter due to the partial cooling. 

The center panel shows the effects of full heating (i?b — 
i?dm) and blowout, as if there is no cooling. For the two kinds 
of galaxies, the baryonic A' distribution becomes signifi- 
cantly higher then that of the dark matter, with Aq = 0.051 
for the bright galaxies and a very high Aq = 0.067 for the 
dwarf galaxies. 

Finally, the right panel shows the result of applying 
the full feedback scheme, considering cooling, heating and 
blowout together. The baryons in the bright galaxies now 
have spins comparable and even slightly higher then their 
dark- matter halos, with Aq = 0.040 (and cj\ = 0.43), while 
the baryons in dwarf galaxies now have spins higher by ~ 
50%, with a;, = 0.068 (and ax = 0.61). 

Thus, the baryonic spin in dwarf galaxies, which are 
made by mergers of small satellites, is dominated by the 
blowout of the gas from these satellites, and it ends up with 
A(, > Ajj^. For bigger galaxies, which are largely made of 
bigger satellites, the dominant effect preventing the over- 
cooling spin catastrophe is the heating, with some contribu- 
tion from blowout. This yields a baryonic spin distribution 
similar to that of the dark matter halos, in general agree- 
ment with observations. 



5 MODEL VERSUS OBSERVATIONS 

We now confront the model predictions with observations, 
in particular those of BBS. We first calibrate the value of the 
one free parameter Vpa by the observed gas deficit in dwarf 
galaxies. We then address the distribution of the baryonic 
spin parameter in these galaxies, as well as the angular mo- 
mentum profiles within dwarf and bright disks. 



5.1 Observations 

Generally, galaxy formation modelers have compared their 
predicted disk sizes to the distribution of di sk scale lengths 
found in large samples o f late type galaxies ( [Courteau 199(: ; 
Mathewson et al. 1992|). A more detailed comparison has 



been performed by ie Jong & Lacey (200C) where they con- 
structed the bivariate space density of galaxies as a function 
of luminosity and size. This allows the spread in observed 
disk sizes to be compared to models and they found that 
the spread corresponded to a a\ of only 0.36 ± 0.03. Ideally 
one would like to compare the actual angular momentum 
in models with observed galaxies instead of just properties 
derived from the angular momentum. 

The only observational attempt so far to measure spin 
parameters for galactic disks is by BBS, which therefore 
serves as the main target for our modeling effort in this 
paper. BBS used fits to the rotation curves and assumed 



an NFW profile for the halos (Navarro et al. 1995) in or- 
der to determine the halo virial quantities J?vir and Vvir for 
14 dwarf galaxies, with an average Vvir — 60kms~^. They 
also assumed a mass-to-light ratio of unity in the R band, 
but their results are not very sensitive to this assumption 
because dwarf galaxies are dominated by their dark matter 
halos. BBS then determined for each galaxy the baryonic 
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Figure 7. The probability distribution of the fraction of mass in 
baryons /d for the data of BBS (shaded), in comparison with our 
model predictions with Vfb = 95kms~^ for dwarf galaxies (left) 
and for bright galaxies (right). The model prediction for the dwarf 
galaxies, with significant blowout, is in good agreement with the 
BBS data. The bright galaxies retain most of their baryons. 




0.05 0.10 0.15 0.20 

Figure 8. The distribution of A' in the sample of dwarf galaxies 
by BBS (shaded histogram) compared with our model predictions 
with Vfb = 95kms~^ for the baryons in dwarf galaxies (curve on 
the right), and the simulation result for dark halos by BD (curve 
on the left). The inclusion of blowout produces a A' distribution 
in good agreement with the data. 



spin parameter Ab which we have converted to Aj, using 
their values for Cvir and equation ^. This leads to an av- 
erage value of Ab ~ 0.063, significantly larger than that of 
simulated dark-matter halos. They also estimated the ra- 
tio fd of disk mass (stars -|- gas) to dark-matter mass, and 
found an average of fd — 0.04, abou t a factor of 3 sm aller 
than the universal baryonic fraction ( Tytler et al. 1999 ) . Fi- 



nally, they measured the j distribution in each galaxy, and 
confirmed the mismatch of angular-momentum profiles. We 
discard one galaxy (UGC4499) from the sample because its 
rotation-curve fit yields strikingly anomalous results. 



5.2 Baryon Fraction 

The distribution of fd values for the dwarf galaxies of BBS 
is displayed in Fig. M as a shaded histogram. These values 
are significantly lower than the generally adopted univer- 
sal value of /b — 0.13, which suggests significant baryonic 
mass loss from these objects. Shown for comparison are the 
corresponding model predictions for dwarf and bright galax- 
ies. We have chosen a value of Vfb — 95kms~'^ in order to 
match the mean of the BBS measurements for dwarf galax- 
ies, fd — 0.04, but one can see in Fig. |^ that the resulting 
spread in fd values is also in good agreement with the dwarf 
galaxies data. For bright galaxies, fd is typically lower than 
the universal value by only ~ 30—40%, refiecting the limited 
fraction of small merging satellites who lost their gas.^ It is 



^ There is one dwarf galaxy with an fd value much higher then 
the rest of the distribution. This galaxy has the lowest virial ve- 
locity in the sample, V^ir = 43kms^^, and the second highest 
concentration parameter, Cvir = 31, suggesting that there may be 
something unusual about this galaxy, either in its history or in 
the BBS fit to it rotation curve. 



encouraging that the obtained value for Vfb is in the range 
of expected values for supernova- driven winds (DS). 



5.3 Spin Parameter of Baryons 

Next, we compare predicted and observed distributions of 
spin parameters for dwarf galaxies. The distribution of Xi 
for the BBS disks of dwarf galaxies is plotted in Figure H 
as a shaded histogram. Shown first in comparison is the 
A' distribution for dark halos in cosmological ACDM sim- 
ulations from BD. Contrary to comments made by BBS, 
these two distributions are not in agreement — the ob- 
served spins are significantly higher, with an average of 0.063 
compared to 0.035 in the dark-matter simulations. Finally 
shown in Figure ^ is our model prediction for the baryonic 
spin distribution in dwarf galaxies (Vvir = 60kms~^), with 
Vfb = 95kms~^. The effect of blowout brings the baryonic 
spin distribution into good agreement with the observed 
sample. Note that this distribution is quite different then 
that of the dark matter, with Aq — 0.068 and a\ = 0.61. 
We find it remarkable that the agreement with the observed 
spin distribution follows automatically from the adjustment 
of V{b to match the observed baryon fraction in dwarf galax- 
ies. In the next section we demonstrate that the result is also 
independent of the details of the feedback model used. 

In the context of the scatter issue mentioned above we 
note that the predicted spread in A' values for baryons in 
bright galaxies is smaller than the spread for the dark-matter 
halos (cta = 0.43 versus 0.5), and is more like the spread ob- 
served by ie Jong & Lacey (2000). We note also that the 



high baryonic spin values obtained in our model for both 
dwarf and bright disks suggests that one does not need to 
eliminate low-spin halos in order to ma ,tch the observed dis- 
tribution of disk sizes [e.g., as done in Mo et al. (199^)]. 
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Figure 9. Examples of baryonic specific angular-momentum profiles for the BBS data (left column) and for our model dwarf and bright 
galaxies (center and right columns respectively). Also drawn is the average (/i = 1.25) dark-matter jr-profile seen in the simulations of 
BD (dashed curve). The BBS profiles clearly lack the high and low j tails seen in the simulations. Likewise, many of our model dwarf 
galaxies also lack this material because of gas loss. The bright galaxies, which retain most of their baryons, have j profiles more similar 
to the dark matter. However, in some cases they are consistent with the profile of an exponential disk with a flat rotation curve (solid 
curve in right column). 



5.4 Angular-Momentum Profile of Baryons 

BBS also measured the specific angular momentum profiles 
of their dwarf galaxies, which they find to be quite different 
from the profiles of dark halos found by BD. The dwarf 
galaxies in their sample show a systematic behavior, with 
a low baryonic fraction, and a significant deficit of angular 
momentum at the two ends of the distribution compared to 
the halo j profiles of BD. Examples of some of the j profiles 
from their data are shown in the left column of Fig. ^. 

We construct a baryonic j profile in each of our model 
realizations, following the sam e method used to produce 
dark-matter j profiles in §2.2, but now including feedback 



and gas loss. Examples of such model profiles for dwarf 
galaxies of Vvir ~ 60kms~^ are shown in the central col- 
umn of Fig. ^ The low-j tail of the halo distribution is 
missing in the baryons in many cases, as expected, because 



the blowout preferentially removes gas from small satellites 
— those satellites that merge to form the low-j halo ma- 
terial. The high-j tail tends to be reduced in the baryons, 
because this tail is often the result of a small satellite that 
comes in with its orbital angular momentum aligned with 
the halo spin, and now has lost its gas. While there is gen- 
eral agreement between the model dwarf galaxies and the 
BBS data, sometimes the model galaxies show a spike of 
low j material that is not seen in the data. We can think 
of several alternative reasons for this apparent discrepancy. 
Either the predicted spike corresponds to a baryonic com- 
ponent that BBS fail to observe (such as faint halo stars), or 
these galaxies which show a spike do not become disk dwarf 
galaxies and are therefore missing from the BBS sample, or 
our model is missing a certain element that should eliminate 
the very low-j spike more efficiently. 

The right column of Fig. ^ shows the model predictions 
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Figure 10. Probability distribution of the spin-profile shape pa- 
rameter f = jtot /imax (the BBS alternative to fi of BD) for the 
dwarf-galaxy data of BBS (shaded histogram) in comparison with 
our model predictions for the baryons in dwarf (thick histogram) 
and bright (thin histogram) galaxies. Shown also is the result for 
simulated halos by BD (smooth curve). The model prediction for 
baryons is in reasonable agreement with the data, though the pre- 
dicted spread is somewhat larger. Bright galaxies show a slightly 
shifted distribution to higher values of (. 

for a sample of bright galaxies, Vvir ~ 220kms~^. They re- 
tain most of their baryons, so their profiles are less affected 
by blowout, and they are more similar to the dark- matter 
j-profiles seen in the simulations of BD. However, in some 
cases they show a reduction of both low and high j material 
for the same reasons that it is seen in the dwarf galaxies. 
For comparison, we also show in the right column of Fig. ^ 
the j profile of an exponential disk with a fiat rotation curve 
truncated at 4.5 disk scale lengths. The j-profile mismatch 
problem in bright galaxies is partly remedied by the com- 
mon presence of low-j bulge components in such galaxies. 
The fact that sometimes the low-j material is lost even from 
bright galaxies allows some of them to become pure expo- 
nential disks, as observed, without the need to hide any low-j 
material in a bulge. 

Fig. |lo| shows the distributions of the quantity which 
characterizes the shape of the j-profile. It compares the ob- 
served values in the BBS dataset with our model predictions 
for the baryons in dwarf and in bright galaxies. Given the 
crudeness of our model, the agreement between our model 
for dwarf galaxies and the data is quite reasonable. The 
somewhat smaller spread in the data compared to the model 
predictions may suggest that either selection effects or gas 
processes tend to favor a particular value of ^. Also shown 
for comparison is the distribution of ( for dark halos as mea- 
sured by BD in simulations. One sees that the values of 
for bright galaxies are slightly shifted to higher values then 
the dark halos. For dwarf galaxies the shift is a significant 
effect, refiecting the mismatch between the j-profiles of dark 
halos and observed galaxies. 



Table 1. The dependence of our results on the model parameters. 
For a given 72 the scale Vf-b is set by requiring the average 
in dwarf galaxies (V^ir = 60kms~'^) to be 0.4. Then we find 
the lowest value of 71 that creates a spin distribution in bright 
galaxies (Vvir = 220kms~^) with a median > 0.035. 
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1.5 


0.4 


0.067 


0.61 


0.9 
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0.64 
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0.036 0.42 



6 ROBUSTNESS OF THE FEEDBACK MODEL 

In this section we examine the sensitivity of our results in the 
previous section to our extremely simple model of feedback. 
We explore a range of values for the exponents 71 and 72. 
For each choice of exponents, we determine Vtb such that the 
mean value of fd for the dwarf galaxies is the same as the 
BBS value of 0.04 We find that the results of § | for dwarf 
galaxies remain practically unchanged. The results for bright 
galaxies do depend on the value of 71, though a value can 
always be chosen such that the median spin for the bright 
galaxies is at least as high as that of the dark halos. 

In Table |l| we show a range of values of 72 and the corre- 
sponding values of Vfb necessary to give our dwarf halos the 
right baryon fraction fd- Then we list the largest value of 71 
for which bright galaxies have a median A[, > 0.035, as sug- 
gested observationally. A maximum value of 3.0 is considered 
for both exponents as higher values produce indistinguish- 
able results. The resulting values of fd, X'o and a\ are listed 
for the two types of galaxies. We indeed see in the table that 
varying the exponents 71 and 72 has a small effect on the 
resulting baryonic spin distribution for dwarf galaxies. We 
conclude that any feedback recipe in which an appropriate 
amount of gas is being ejected from dwarf galaxies should 
resolve the angular-momentum crisis for dwarf galaxies in 
a way similar to our results in The angular momentum 
catastrophe for bright galaxies is prevented if 71 is suffi- 
ciently small. 

If 72 is as small as 0.5, as in the first row of Table |l|, 
then the required value of Vfb needed to fit the mean value of 
fd in dwarf galaxies is very high (185kms~^). This velocity 
is significantly higher than the crude upper limit allowed by 
the energetics (DS). Also, for such a large Vfb value bright 
galaxies have a median A' much higher then the dark matter 
even for the maximum value of 71. Thus, the model yields 
sensible results only if 72 is limited to the range 0.8 — 3.0. 

The shaded region in Fig. ^ shows the range of curves 
corresponding to the range spanned in Table |l| We note that 
a tight constraint is enforced by the data on this range of 
feedback models, that halos of ~ 35kms~^ must loose about 
half of their baryons. 



7 CONCLUSION 

We used a simple model to address the two issues that cre- 
ate the angular-momentum crisis of galaxy formation within 
the CDM scenario. One, the angular momentum catastro- 
phe is that the spin of the baryonic component in galax- 
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ies is typically comparable to that of the dark-matter halo 
(in bright galaxies) or even larger (in dwarf galaxies), while 
simple theoretical arguments and simulations involving gas 
dynamics predict a significant spin loss by the baryons due 
to over-cooling in merging halos. The other, the mismatch 
of angular-momentum profiles, is that the baryons in each 
galaxy tend to lack the low-ji tail (and the high-ji tail) of the 
distribution as predicte d by simulations fo r the dark halos. 

In an earlier paper (Mailer et al. 2002), we showed that 
a simple algorithm, based on adding up the orbital angular 
momenta of the mergers in random realizations of merger 
historie s, can successfully reproduce the distribution of spins 



to incorporate a more sophisticated feedback recipe into the 
model using the full machinery of semi-analytic models of 
galaxy formation. This will be a step towards the long-term 
goal of implementing feedback in full-scale hydrodynamical 
cosmological simulations. 

Our work is an attempt to resolve the angular- 
momentum problems within the standard framework of 
CDM cosmology, which is so successful on large scales, using 
the effects of feedback which we know must occur. Another 
approach is to appeal to a different cosmological scenario, 
WDM, wher e the dark-matter particles are "warm" rather 
than "cold" ([Hogan 1999| ; [Hogan fc Dalcanton 200o| ; [PageL 



among dark- matter halos as measured in N-body Simula- fc Primack 1982| ) . WDM is less robust than CDM because it 



tions of the ACDM cosmology. We showed here that an ex- 
tension of such a model also reproduces the characteristic 
angular-momentum profile, i.e., the distribution of specific 
angular momentum within halos. This provided the basic 
tool for addressing a possible resolution to the spin prob- 
lems by incorporating the effects of supernova feedback on 
the gas in halos before they merge into bigger halos — a 
process we have termed spin segregation. 

A simple analysis of how the orbital angular momen- 
tum in a merger turns into a spin profile suggested how 
feedback effects in the satellite before the merger event can 
eliminate the problem of spin loss. The effective size of the 
gas component within the incoming satellite determines its 
tidal stripping position in the halo and thus the final spin 
that it will be left with after the merger. The finding, us- 
ing the orbital-merger model, that the low-end tail of the j 
distribution originates in many minor mergers, that tend to 
cancel each other's angular momentum, provided a possible 
solutior i to the spin-profile mismatch problem. The blowout 



requires fine-tuning of a new parameter, the particle mass, 
which is constrained to be ~ 1 keV. Still, it is worth inves- 
tigating since it may remedy the angular-momentum prob- 
lems without appealing to strong feedback effects. The main 
distinguishing feature of the WDM scenario is that the for- 
mation of small halos is significantly suppressed, such that 
the validity of the explicit picture of halo buildup by the hi- 
erarchical build up of small halos becomes limited. Despite 
t his difference from CDM, an N-body simulation of WDM 
(Bullock et al. 2001) indicates that the angular-momentum 



properties of halos remai n basically unchang ed. This is not 
very surprising because in Mailer et al. (2002D we found that 
the same angular-momentum properties can also be inter- 
preted as a general result of tidal-torque theory, indepen- 
dent of the explicit picture of mergers (see also BD). In 
the absence of small halos at early times, one may expect 
less over-cooling in halos before they merge into other halos 
and thus less angular-momentum loss by the ba ryons. In- 



deed, hydrodynamical simulations of this scenario (Sommer 



of gas f |:om small incoming halos would eliminate the low-j Larsen fc Dolgov 200 1| ) indicate that the angular momentum 



tail of the baryon distribution in the merger product, as ob- 
served. The blowout has a particular strong effect in dwarf 
galaxies, because they are made of smaller progenitors that 
tend to loose more of their gas. This results in a higher 
spin parameter for the baryons than the dark matter, as 
observed. 

We constructed a simple semi-analytic model for simu- 
lating this process using a simplified model for the effects of 
feedback as a function of halo mass, including heating and 
blowout. For a given choice of exponent 72, the model has 
one free parameter, the characteristic halo virial velocity Vfh 
for which the energy inputted to the IGM by supernova is 
sufficient to counter the effect of cooling. By matching the 
low baryonic fraction in the dwarf galaxies observed by BBS, 
we found that for 0.8 < 72 < 3 the characteristic velocity has 
to be in the range 60 < Vfb < 130kms~^, which falls within 
the range of theoretical predictions (e.g., DS). We then found 
considerable agreement between the model predictions and 
the observed data for both the distribution of the spin pa- 
rameter and the angular-momentum profile of baryons in 
dwarf and bright galaxies. We also noted that the same ba- 
sic model may explain the other unresolved issues concern- 
ing angular momentum in galaxies, such as the spread in 
observed disk sizes and the identification of halos that form 
late type galaxies. The success of the model in matching 
several independent observations indicates that this simple 
model indeed captures the main features necessary for a full 
scenario involving feedback and mergers that can resolve the 
spin crisis in more detail. The next natural step should be 



catastrophe is significantly remedied. However, the angular- 
momentum profile mismatch is still expected to be valid in 
WDM, and, in the absence of small halos, the feedback ef- 
fects are weaker and may not be enough for resolving the 
problem. Thus the j-profile mismatch may be a crucial dis- 
criminator between solutions to the problems of CDM. 

It is clear that some sort of heating may provide the cure 
f or another problem of CDM, the dw arf satellite problem 
( Klypin et al. 1999| ; Moore et al. 1999), where the predicted 
number of dwarf halos is much larger than the observed 
number of dwarf galaxies. 



Bullock et al. 
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demon- 
strated that cosmological photo-ionization due to feedback 
from UV sources such as early stars and quasars can solve 
this problem. Feedback from supernova would qualitatively 
have a similar effect. While the number of dwarf satellites is 
automatically suppressed in WDM, it seems that the inclu- 
sion of minimum realistic feedback effects would reduce the 
predicted number of galaxies to significantly below the ob- 
served number, and thus be an overkill (J. Bullock, private 
communication) . 

Furthermore, it is becoming clear (Dekel et al., in prepa- 
ration) that the key elements of our model, namely the tidal 
effects in mergers and the feedback effects in small halos, are 
also very relevant in understanding and resolving the third 
problem of CDM. This is the cusp/core problem, where the 
halos in sim ulations typically sh o w steep cusps in t heir in- 
ner profiles ( Navarro et al. 199£| ; Moore et al. 199£| ), while 
observations indicate ffat cores at least in some low-surface- 



brightness galaxies ( de Blok et al. 2001 ) . An analysis of tidal 
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effects explains the necessary formation of an asymptotic 
cusp in halos as long as satellites continue penetrating into 
the halo center. Feedback effects may puff up small satel- 
lites and prevent them from penetrating the core, thus al- 
lowing a stable core. These studies together indicate that 
our model indeed grasps the relevant elements of the com- 
plex processes involved, and that feedback effects may in- 
deed provide the cure to all three major problems of galaxy 
formation in CDM. 
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